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Keywords: Vector meson, diffraction, QCD. 



1. Introduction 

Diffractive vector meson photoproduction at large momentum transfer, 7p — ► VX , is an 
interesting process experimentally, because of the clean signal: an isolated vector meson 
with large transverse momentum, that is separated from the proton remnant by a large 
rapidity gap. Theoretically, it is advantageous because the large momentum transfer \t\ 
provides the hard scale necessary for a perturbative QCD description of the process in 
terms of hard colour singlet exchange. It has therefore been proposed as an ideal testing 
ground for BFKL dynamics |l[ |2| . 

This process was recently investigated at HERA pi |J. The ZEUS analysis includes 
the angular distribution of the decay products, giving access to helicity amplitudes of the 
photon-meson transition. The measurements gave somewhat surprising results. Both for 
p and 4>, the dependence of the cross-section on the momentum transfer t is approximately 
power-like, \t\~ n , with the exponent n ~ 3. It can be deduced from the data that the 
leading contribution comes either from the process with no helicity flip or from the one 
with double helicity flip. The subleading helicity amplitudes were shown to be one order 
of magnitude smaller than the dominant one. 

An analysis based on leading-order perturbative QCD predicts a steeper decrease of 
the cross-section, n ~ 4, for the no-flip and double-flip amplitudes ||. The single-flip 
amplitude M + o gives a differential cross section with n ~ 3, and should become dominant 
at very large momentum transfer \t\. It is clear that these properties are in contradiction 
with the experimental results. 

Improved theoretical understanding came about after the suggestion that the real pho- 
ton may couple to a quark-antiquark pair with a significant chiral odd component ||. In 
general, the quark-antiquark pair can couple with either chiral even or chiral odd compo- 
nents, the latter vanishing in perturbation theory in the limit of massless quarks. However, 
it is not clear that the mass should be interpreted as the current quark mass. The answer 
to this question turns out to be essential for the phenomenology. Estimates based on QCD 
sum rules were performed in giving a characteristic value of the non-perturbative mass 
parameter rather close to the constituent quark mass for light vector mesons. At moderate 
\t\ ~ 10 GeV 2 , the chiral odd amplitude with no helicity flip (M°| d ) is then expected to 
dominate. In this paper we compute the chiral odd and chiral even amplitudes, leaving 
our results explicitly in terms of the quark mass m. An alternative explanation for the 
dominance of M ++ and the emergence of dimensional scaling can be found in ||. We shall 
return to discuss the connection of our approach with that of Q in our second paper, where 
we shall focus on the phenomenology Q. 

The analysis of ||] was restricted to the lowest order approximation for colour singlet 
exchange, that is the exchange of two gluons. In fact, the diffractive production of vector 
mesons occurs at energies y/s much larger than the momentum transfer. In such kinematics, 
the perturbative QCD corrections to two gluon exchange are enhanced by powers of large 
logarithms of the energy, log(s/|i|). The leading logarithmic corrections are proportional 
to [a s log(s/|i|)] n at the nth order of the perturbative expansion. Thus, it is not sufficient 
to consider the lowest order approximation. A convenient way to perform the resummation 
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Figure 1: Vector meson photoproduction at large momentum transfer via hard colour singlet 
exchange. 



of leading logarithmic terms to all orders in a s is given by the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) equation ||. This approach turned out to be successful in the case of 
diffractive production of charmonia at high \t\ [Q, [2|, ||, [ll]]. Interestingly enough, the 
t-dependence of vector meson production was successfully described in |l(J by a BFKL fit 
with a non-relativistic wave function for both light and heavy mesons, we shall comment 
further on this finding later. 

The purpose of this paper is to present analytic results for the chiral-even and chiral- 
odd helicity amplitudes with complete leading logarithmic BFKL resummation. This cal- 
culation will give control of important QCD effects in vector meson photoproduction at 
high energies. A detailed comparison to available data will be performed in a forthcoming 
paper @. 

The paper is constructed as follows: The basic ideas are introduced in Section 2, the 
helicity amplitudes at the lowest order are given in Section 3 and their BFKL evolution is 
found in Section 4. Brief conclusions can be found in Section 5. 



2. Hard colour singlet exchange 

Let us consider a frame in which the photon-proton collision occurs along the z axis. The 
momentum transfer vector q is dominated by its transverse part q, t = q 2 ~ — q 2 . The 
possible helicity states of the incoming quasi-real photon are characterized by transverse 
polarization vectors 

e ± = T^(l,±i), (2.1) 

and similarly for the transversely polarized vector mesons. For the meson, though, the 
longitudinal polarization e° is also allowed. Thus, the possible photon-meson transitions 
may be described using three independent helicity amplitudes M ++ (no-flip), M + q (single- 
flip) and M+_ (double-flip) 1 . 

1 There are also the corresponding amplitudes M , M_o and M_+ which satisfy M++ = A/__, M-\ = 

M_+ and M +0 = -M-o- 
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The diffractive process jp — > VX at large momentum transfer t (see Fig. |l]) takes place 
by exchange of the BFKL pomeron. It has been demonstrated that at large momentum 



transfer, the hard pomeron couples predominantly to individual partons in the proton [12]. 
Thus, the cross-section may be factorized into a product of the parton level cross-section 
and the parton distribution functions, 



da{jp — > VX) ( 42V, 



G(x 3 ,t) + Y}q f (^t)+q f (x 3 ^ Miq ^ Vq \ (2-2) 



dt dxj V ( N c ~ 1) 

where N c = 3, G(xj,t) and qf(xj,t) are the gluon and quark distribution functions respec- 
tively, and W 2 is the jp centre-of-mass energy squared. The struck parton in the proton, 
that initiates a jet in the proton hemisphere, carries the fraction Xj of the longitudinal 
momentum of the incoming proton. For future use, we introduce the integrated quantity 

/da( 7 q^Vq)\ ! MnP^l. (2 . 3) 

\ dt I J cuts dtdxj 

The partonic cross-section, characterized by the invariant collision energy squared 
s = XjW 2 is expressed in terms of the amplitudes M\\r(s,t), 

^ = ±- {\M ++ (s,t)\ 2 + \M +0 (s,t)\ 2 + \M + .(s,t)\ 2 ) . (2.4) 

In the case of an unpolarized lepton (photon) beam, it is impossible to measure directly 
the helicity amplitudes due to interference effects. Thus, in the experiment one selects a 
sample of vector mesons and measures the angular decay distribution of the decay products 
H. The following spin density matrix elements can be determined: 

r °°" (|M ++ |2 + |M +0 | 2 + |M + _|2)' W 



04 1 (m ++ m; + m + -M* ) 

7-10 2(|M ++ | 2 + |M +0 | 2 + |M 4 l 2 \' ' 



04 _ (TZe(M ++ M* + _)) 



- n dM+^ + lM+of + lM+.l 2 )' 

where (...) denotes the integration of the parton level quantities over partonic Xj with the 
appropriate cuts (c.f. (fO|)). 



3. The lowest order amplitude 

In the high energy limit it is convenient to use the dipole representation for the calculation 
of the scattering amplitude. In this limit the helicities of quarks, the large light-cone 
components of their momenta, and the transverse positions, are conserved in the scattering 
process. The small components of light-cone momenta are integrated out. Thus, the impact 
factor corresponding to a transition of a photon with helicity A to a meson with helicity 



-3- 



A', due to the coupling of two gluons with momenta k and q — k, may be represented in 
the following way 



$7(A)V(A') (fe) g) = j duj d 2 r [** v(y) (u,r)UT(k,q;r,u) [^ {x) ( U ,r)) af3 , (3.1) 

where Vl/y (^-y) denotes the vector meson (photon) wave function, u (u = 1 — u) is the 
quark (antiquark) longitudinal momentum fraction, and r the transverse quark-antiquark 
separation vector. Summation over spinor indices a, and (3 is performed. T(k,q;r,u) is 
the hard amplitude of the QCD dipole scattering, which, up to a constant factor, equals 

jdipole _ e iqru (\_ e ^ ikr \ (\ — e - i (l^ k ) r \ ; (3.2) 



see H]. If we consider the meson photoproduction off a quark, the complete amplitude 
takes the form: 

M X x(q) = J d 2 fc fc2 <F^ v ^(fc,<?) <& gg (fc,g), (3.3) 

where the quark impact factor Q qq (k,q) is a constant. We choose the normalization 2 

$ 9 "(fc,g) = l. (3.4) 

In order to produce the meson, we follow [|j in assuming that the quark-antiquark pair 
has small transverse size (as determined by the exchanged transverse momentum which is 
shared approximately equally between the two gluons (l|). This approximation allows us 
to expand the relevant hadronic matrix element about the light-cone, keeping the leading 
terms in x 2 (where x is the space-time separation of the quark and antiquark) . The relevant 
light-cone wavefunctions are then taken from the QCD analyses of |l3|, 14, H|. 

The photon wavefunctions in the momentum representation were obtained using the 
perturbative approach in ||16|| . These wavefunctions contain spinorial structures with defi- 
nite chiral parity. The chiral even components (corresponding to vector 7^ and axial vector 
7^7 5 coupling to quarks) yield non-zero contributions in the massless limit, and the chiral 
odd ones (tensor cr^ u ) are governed by the quark mass m, and vanish when m — > 0. Thus, 
it is convenient to classify the contributions to the impact factors and amplitudes according 
to their chiral parity. Then we have 

*7(A)V(A') = ^(A)V(A') + ^(A)V(Y) ^ 

and 

M\\i = + M$f, (3.6) 

for all photon and meson helicities A = +, — and A' = +, — , 0. In || the available informa- 
tion on the vector meson wave functions was employed to derive both the chiral even and 



2 Only the normalization of the full amplitude is relevant in our calculation, so we move all non-trivial 
factors into the photon-meson impact factor. Caution should therefore be exercised if these impact factors 
are to be used elsewhere. 
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the chiral odd components of the amplitudes (p . 3|) . These amplitudes, given at the leading 
order, are the starting point for our BFKL analysis. 

The chiral even impact factors for 7(A)V(A') computed in ||] are given by 

^ e + „ )y(0) = ^|^^m^ 1 Hr|)^^/ dipole /y(l-2 U )0 || (n), (3.7) 
= c v J m\r\ Kx(m\r\) (e« • e (+)*)/<» pole 

xfvMy^ (]± hiv)+ j% Mv ]\, (3.8) 

\0 u / 

•jffi V <-»- -Or /•j^ mglH r|) (6 ' +Kr ; ( f" ■ r) /^« 
J 47T \r\ 

xf v M v (u 2 j d ^^(v)+u 2 f (3.9) 

V u I 

with 

iV 2 - 1 

Cv = io? s ^-eQ v , (3.10) 

where m is the quark mass, fy is the meson decay constant, and we take the asymptotic 
wave function (f>\\(u) = 6u(l — u) for the light vector mesons. The effective quark charges 
in the vector mesons are given by Q p = l/y/2, Q u = 1/3V2, and Q = -1/3. 

The chiral odd impact factors for 7(A)V(A') with a perturbative (QED) photon wave- 
function read 

^MV(0) = _ lCy J fL^ mKQ{mlrl) (6 (+) . r) /dipole 

xtfMvuuU ^<p ± (v)- j ^MvU , (3.11) 
*2S )V(+) = C,|^ ra ^Hr|)( £ W. e W*)/^/^i(n) 1 (3.12) 

^ )VH = ^cv/^mtfoHrl)/^ 

x ( e (+) • r)(e(-> • r) My 6(uu) 2 , (3.13) 

where <t)±(u) = 6uu, and the tensor decay constant fy are introduced, and we used the 
asymptotic form of distribution functions h\\(u), hs(u) and 4>±(u) in ( |3.13 ) [||, 13]. 

In what follows, the complex representation of the two-dimensional transverse vectors 
will be used, e.g. r = r x + ir y . This leads to the substitutions reW = —r/y/2, d 2 r = d 2 r = 
drdr* /2i, |r| — > |r| etc. 

Thus, at the lowest order, the chiral even amplitudes for vector meson photoproduction 
are given by 

M +° en= * Cv 7 i?i^ d ^ m v\ KiMr ^ (3 - 14) 
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M:T = Cl 7 Hr|*.(m|r|)/**A*v 



\0 



d 2 fc d 2 r r 2 
fc-q) 2 ~1 

%,|W+n 2 / , (3.16) 



x I it 2 ' u " ^ ' - 2 



u 

and the chiral odd ones 

M$ d = ^ f ^Jl^L fL^ mrKo (m\r\) / di P° lc 

xV2tfM v uu (j j^<t> x (v) \ , (3.17) 

\u / 

Mf + d = CV / fc2( f fc g)2 ^^m^ (m|r|)/ di P° lc tf^u), (3.18) 

M odd = j *k fL^ mr 2 KQ{mlrl)fiP 0le f T M 2 (3 ^ 

o J A; (fe — g)^ 47r 

Let us consider the behaviour of the amplitudes in the massless limit. For all chi- 
ral even amplitudes, one obtains non-zero results in this limit, as mKi(m\r\) — ► l/\r\ 
for m — > 0. The situation is different for the chiral odd amplitudes. The factor of 
mKo(m\r\) ~ — mlog(m|r|) implies the vanishing of all perturbative chiral odd ampli- 
tudes in the massless limit. However, as discussed in [||, chiral symmetry breaking effects 
can induce a sizeable chiral odd coupling of the photon to quarks. In particular, the 
characteristic mass parameter of the chiral odd coupling was estimated in || to be about 
2-K 2 fj/N c ~ 0.46 GeV. Here we account for chiral symmetry breaking by keeping in mind 
that the quark mass parameter need not be small. We refer to Appendix B for a more 
detailed discussion of the limit m — > 0. 

We should also remark that we are able to reproduce the original formulae of ju]] by 
considering the only surviving amplitudes in the limit that the quark and anti-quark are 
collinear, i.e. M° d f and M^f 1 . We also note that if we make the further constraint that the 
quark and anti-quark share the meson's momentum equally, we obtain the formulae of PJ 



for the M_|__|_ amplitude, the only surviving amplitude in this case. In [1C], this amplitude 
was shown to be in good agreement with the ZEUS data (after including BFKL effects) 
allowing us to conclude that a sizeable chiral odd contribution seems to be required by the 
data. 
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4. The exact BFKL amplitude 
4.1 Generalities 

The BFKL kernel in the leading logarithmic approximation exhibits, in the impact pa- 
rameter representation, invariance under conformal transformations [jl7). The conformal 
symmetry of the kernel permits the following expansion of the amplitude in the basis of 



eigenfunctions E n>v [17]: 

M AB (z, q ) = W E j_ 



v 2 + n 2 /4 

dv 



_ x [^ 2 + (n-l) 2 /4][^ 2 + (n+l) 2 /4] 
where 

Xniy) = 4fte(^(l)-^(l/2 + |n|/2 + «/)^ (4.2) 
is proportional to the eigenvalues of the BFKL kernel and 



2tt V A2 , 

(A is a characteristic mass scale related to My and \t\). 
f d 2 k f 

In,M) = J (2^)2* ( k ' q ^ J d2 Pi d2 p2 E nAPi>P2) e Mik- Pi+i(q-k)- p 2 ), (4.4) 

and analogously for the index B. The eigenfunctions are given by 

\ -2+1/2 / / \ *\ -u+1/2 

Pi ~ Pi \ * ' I ( Pi ~ P2 X X 



V / VV PlP2 ) ) 

where p = n/2 — iv and \x = —n/2 — iv. Here k and q are transverse two dimensional 
momentum vectors, and pi and p 2 are position space vectors in the standard complex 
representation. The functions $> A = Q'YWVix ) an( j _ ^gq are ^ e j m p ac ^ factors. 

The quark impact factor in representation ( |4.4[) was found in generalizing the 



Mueller- Tang subtraction [19] to non-zero conformal spin; 

M W r(l/2 + n/2 + ii/) l4 '°' 

for even and 7^ = for odd n. 

We now compute the impact factors for the A = j(X)V(\') transition for all helicity 
states and all conformal spins n. The following integrals are to be evaluated 



/d 2 k a f 
J^y$ ' (k,q) J d 2 p l d 2 p 2 E ntU (p 1 ,p 2 ) 



x exp(ife*pi/2 + ikpl/2 + i(q* - k*)p 2 /2 + i(q - k)p* 2 /2). (4.7) 
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Changing variables to p\ = R + p/2, p2 = R — p/2 and integrating over d 2 k we get 

tiM = J d 2 p^( P ) J d 2 RE ntU (R, P yi R , (4.8) 
where &q(p) is the Fourier transform of the impact factor in the momentum space <& A (k, q), 

*tip) = J -0^^ A (k,q)e^-^. (4.9) 



The integration over R in Q4.8| ) has been done by Navelet and Peschanski pOf| , resulting in 

lL(Q) = J d 2 p^(p)E^ u (p), (4.10) 
where the mixed representation eigenfunction En,v(p) is given by 

= h n-y Of) ^ (7) W2 r ^ " ^ " n ^ T{l ~ iv + n / 2 ) 

x |p| [J^p/A) J~,(qp*/4) - (-IT J-,>(q*p/4) J-~M*m (4.11) 



and 

2^vr 3 r(|n|/2-^ + l/2)r(|n|/2 + ^) 
n,v |n/2| -«/ r(|n|/2 + ii/ + l/2)r(|n|/2-ti/)" l ' ' 

The integral over fc of impact factors (^^, |3.7| , 3.9 ) and (|3.12 , 3.11] , |3.13[) may be easily 
performed since k appears only in y dl P ole . Due to the resulting 5 functions S(p), 5(r — p) 
and 6(r + p) the integration over r is trivial. Using the fact that En,u{p) = for p = 0, and 

KM = (-V n E%(-p) (4-13) 
one arrives, for all impact factors, at complex integrals of the type 

I a p(v,n,q,u;a) = mj d 2 p p a+1 p*? +1 K a (m\p\)e^ 

x [J^q*p/A)J fi (qp*/A) - (-1)™ J^(q* p/4)J^(qp* (4.14) 

where £ = 2u — 1 and K a {x) is the modified Bessel function, displaying effects of the quark 
mass. The parameter a equals 1 for the chiral-even and for the chiral-odd impact factors 

of §. 

4.2 Amplitudes 

The helicity amplitudes can then be expressed in the following way 



M even = f' ^ ^ _ _ ^ 

4V2 \q\ Jo 



AV2\q\ 

n=+oo oo U 2 + n 2 exp[ X 2n(^)z] , , 

2^ / rft/ T 2 i 7 1 /o\2u 2 i 7 — i 1 /o\2i — S — h-i[y,2n,q : u\l), 

' + (ra — l/2)^J[z/ z + (n + 1/2) J sm(z7rz^) 



(4.15) 
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M evcn = v , du6u ( 1 _ u )( n i_ n+1 / 2 ) 

8 M jo 

n=_^oo Z. 2 + n 2 eMX2n(v)z\ 



n=— oo 



[z/ 2 + (n - 1/2) 2 ] [z/ 2 + (n + 1 /2) 2 ] sin(ivri/) 

(4.16) 



+ m Jo K 1 



n=+oo f00 

x > I du 

-co 



n=— oo 



i/ 2 + n 2 



exp[x2»(^)z] r / 9 n 

[ I /2 + (n-l/2) 2 ][i/ 2 + (n + l/2) 2 ] sin(i7rzy) +1_U ' h 



(4.17) 



where we performed integration over v in (3.14, 3.15| , 3.16). 
Analogously, one may express the chiral odd amplitudes 



Ml 



odd _iC v ftM v f 1 



8y/2\q\ 



du 6 u(l — u) (1 — 2u) 



E 



dz/ 



z^ 2 + n 2 



[z/ 2 + (n - 1/2) 2 ] [v 2 + (n + 1/2) 2 ] sin(mv 



ex-V{X2n{y)z\ , . 
r-r: — : — J I _ i (^, 2n, g, u; 0), 



(4.18) 



Af 



odd _ Cvfy 



du 6 u(l — u) 



E 



(/;/ 



z^ 2 + n 2 



exp[x2n(^)^] Tin n\ 

[„2 + (n- l/2) 2 ][z, 2 + ( w +l/2) 2 ] shuW) J -i-l^.^9»«. u ^ 



(4.19) 



ii^-odd CvfyMy f 1 2 2 

M + _ =^T7-y o **6« (l-n) 

f, ^ 2 + n 2 exp[ X2n (^)z] 

x d V + (n _ 1/2)2][ , 2 + (n + 1/2)2] sin( ^) J i -j 2 ^ ?. u ; o)- 

(4.20) 
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The derivation of integrals ( 4,14| ) for integer difference 3 a — (3 is given in Appendix A, 
and the result reads 

I aP {y, n, q , U ;a) = - ^T(a/2 - C)r(-o/2 - C) r| (i sign (1 - 2u)) a -^ n 

4 \ 4 

1^-1 [sm7r(a + )« + C) B(a, p,q* ,u,() B((3, p,q,u* ,() 

- (-l) n sin7r(a-M + C)-B(a,-//,9*,u,C)B(/3,-?,g,«*,C)] (4.21) 

where we have introduced the dimensionless parameter r g = 4m 2 /|(;| 2 and the conformal 
blocks 



B(a,p,q*,u,0 = (-4nn)-^ +2+a+ ^ 2 fij°2" 



r(/i + 2 + q + C) 

r(/x + i) 



2Fl V 2 ' 2 ^ + 1 ' 4^ J • (422) 

Note, that the sums are performed over even conformal spins 2n due to properties of 
the quark impact factor (^1]) . It turns out that odd conformal spin components of photon- 
meson impact factors Q3.7| , |3.8| , |3.9| ) also vanish. The impact factors are symmetric under 
the coordinate transformation given by (p — * —p, u — > 1 — u). The parity behaviour 
( 4.13 ) of conformal eigenfunctions applied to ( ^4.2l| ) leads to the relation 

I a p(v, n, q, u; a) = (-l) n+a "^ I a p(v, n,q,l- u; a). (4.23) 

These results may be combined to give I%$ V ^ ^ = for odd n. 

The obtained expressions for the impact factors are rather complicated and the deriva- 
tion is rather involved. Therefore we performed an independent test of the results by com- 
paring a direct numerical integration of (4.14) over p to the analytical formula (4.21) for 
some arbitrary sets of parameters, with the quark mass tending to zero. Full agreement 
was found in all cases. 



5. Conclusions and Outlook 

Helicity amplitudes for diffractive vector meson photoproduction off a proton were investi- 
gated using perturbative QCD methods, focusing on scattering occurring at high energies 
and high momentum transfer, for which the target proton dissociates. The process is de- 
scribed in terms of hard colour singlet exchange between a parton from the proton and the 
quark-antiquark component of the photon. 

We re-derived all the relevant impact factors in terms of the meson light-cone wave- 
function. The helicity amplitudes were given in the two gluon exchange approximation 
with arbitrary quark mass and classified according to their chiral parity. We suggest that 
the perturbative formulae with the constituent quark mass reasonably approximates the 
effects of chiral symmetry breaking. 

3 The condition of integer a — (3 is necessary for the amplitudes to be single valued. 
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As the scattering occurs in Regge kinematics, the perturbative corrections are enhanced 
by large logarithms of the energy, and a full resummation of the leading logarithmic cor- 
rections is necessary. Therefore, QCD corrections to the lowest order amplitudes were also 
considered. The amplitudes of the colour singlet exchange may then be described in terms 
of the BFKL equation. Thus, we derived the projections of the photon-meson impact fac- 
tors on the conformal eigenstates of the BFKL kernel for an arbitrary conformal spin. Both 
the chiral even and chiral odd impact factors are obtained in terms of a single complex line 
integral. Using those results, the BFKL evolved helicity amplitudes describing the meson 
photoproduction off a parton were given in an analytical form. These are the main results 
of this paper. 

Our results are correct in the leading logarithmic approximation (LLA). Consequently, 
we are unable to make reliable predictions for the absolute normalization of the cross- 



sections (due to the ambiguity in the scale A). Moreover, next-to-leading logarithms 21, 22 



are well known to be important in the case of the t = BFKL equation (they have yet 
to be investigated in the non- forward case). Attempts to incorporate at least some of the 
next-to-leading logarithms have been performed and indicate that the LLA is probably 
quite reasonable provided one treats the leading eigenvalue of the kernel as an effective 



parameter whose value in LLA is determined by fixing a s ]ll|, |23|]. The fixing of a s is also 



supported by the work of [24, 25]. 

It remains to confront our calculations with the i-dependence of the cross-section and 
the spin density matrix elements which have been measured at HERA ||, |I| for the p, 4> 
and J/ip mesons. This study will be the subject of a second paper Jjj], and should provide 
a stringent test of the BFKL approach, in addition to constraining the photon and meson 
wavefunctions. 
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A. The integral for arbitrary quark mass 



It is convenient to use the variable £ = u — u = 2u — 1 rather than u, and to compute the 
integral for £ an arbitrary complex number and then continue analytically into the physical 
region —1 < £ < 1. We can then write (4.14), using 2i d 2 p = dpdp*, as 

dp dp* 



T± = m 



2i 



P 



a+y /J+l Ka (m\p\) ei«V a +€«*) J±^(q*p/4) J ±Tl {qp*/A) (A.l) 



so thatT = T+ - (-l) n X_. 

Let us start by computing the integral ( A.l ) for massive quarks. The massless case 
will be obtained in a suitable limit. 
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The Bessel functions J are defined with a cut, but the products appearing here are 
single- valued. To make the cut structure of the integral more explicit it is convenient to 
represent the Bessel functions in terms of confluent hypergeometric functions, which are 
individually single- valued, by using the relation p(J 



Ja-i/2(z) = 00, 2a; -2iz) 



Z\0"-l/2 e i 



r((r + i/2) 



(A.2) 



Inserting this into Eq. ( |A.1| ) we have 



64*" q*^ 



r( M + i)r(/i+i) 

x exp 



m 



dp dp* a+1+M » 
2t ^ ^ 



« — : — 9 P + 2 : 9/9 



K a {m\p\) 



4 J ' 4 
X (n + i, 2/i + 1; -%q*p) (# + 5, 2/z + 1; -|gp*) 



(A.3) 



The Bessel function can be replaced with the line integral representation 4 



#a(* ) = - - / — r 

2\X) Jc-ioo 2vrz 



(6/2 + a/2 - z)r(6/2 - a/2 - z) ( — 



x 



(A.4) 



for C < Re[b/2 + a/2], Re[b/2 - a/2]. Substituted into Eq. flOD , we now get the following 
integral: 



m 



T(p + l)r(p + 1) 2 



C'+ioo 



2\ C 



dp dp* Q+ i +M+c * (3+i+jj+c 
. p p e-x-p 



■.e + i * , . r + i * 

* — - A — q p + % — — qp 



x <j> (p + i, 2p + 1; -fg*p) (p + ± 2p + 1; -fop*) 



(A.5) 



where we have changed to the integration variable Q = z — b/2. 

The p and p* dependence of the integrand is clearly factorized, but to disentangle the 
integrals we have to consider carefully the analytic structure. It is convenient to rescale 
the integration variables, 



At 

7 q(C + 1) * 



(A.6) 
(A.7) 



4 This can be seen either by the inverse Mellin transformation or by expressing K a as a Meijer G function 
and using its line integral representation. 
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giving 



qaJu q *n q H m rC'+ioc d( . , /m 2xC 



T(p + l)T(p+l) 2 
dWdZ 



/ ^r(a/2 - c)r(-o/2 - 



w a+l+^+C z l3+l+n+C e -{W+Z) 



2i 



/ 1 2VF \ / 1 2Z \ 

x 0hj + _ 2/x + l;— J ^l/*+2' 2 ^ + 1 ;^piJ • ( A - 8 ) 

Let us now fix Z to be a constant real number. If Z > there is a branch cut along 
the positive real axis from zero to infinity in the complex W plane, and if Z < the cut is 
along the negative real axis. For the integral to converge Z must be positive, and so to get 
a non-zero contribution we have to close the contour around the cut in the right half-plane, 
finding the discontinuity across the cut. Keeping Z fixed, the W integral can be written as 

j> dW{-WZ) u e~ w e~ z F{W) (A.9) 

where u is a non-integer complex exponent and F(W) is the non-problematic, single- valued 
part of the integrand. Let the contour C be the two rays W± = re ±l<5 , where r is integrated 
from to infinity and 5 is a positive real number, anticipating the limit 5 — * 0. We then 
have 

j> dW(-WZ) u e- w e~ z F(W) 

/•oo 

= lim / dr[(-ZW_) u e- w - - (-ZW + ) u e - w +}e- z F(W) 

/•oo 

= 2isin(vrn) Z u e~ z / drr u e~ r F(r). (A.10) 



That is, factorizing the two integrals, we get an extra factor sin(7nt) ensuring that the 
expression is single- valued. We have chosen the direction of the contour that gives the 
correct overall sign. Applying this to ( |A.8| ), with u = a + fi + £, and using the fact that 
fi — Jl = n and a — (3 are both integers, leads to 

U iu q*^q' 1 m f c ' +io ° d( . . . . AN /m 2 \ C 

sin7r(a + // + C) (-i)° HJ+n 



x / 00 dTyiy-+ 1 +^+C e -^ (/) /' /i + I )2/i + 1 .^^ 
Jo V 2 4 + 1/ 

x J™ dZ e~ z <j> (ju + i 2/1 + 1; ^) . (A.ll) 
The last two integrals are instances of the standard integral [ p6( | 

;>oo 

/ due- u u b cf>(a;2a,ku) =T(b + l) 2 F 1 (a,b+l;2a;k), (A.12) 
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and when using this, ( |A.11| ) is reduced to the single line integral 

sin^(Q + ^ + C) (-z) a_/3+n 



/ i 2 
x r(a + 2 + n + C) 2-P1 M + o' a + 2 + A* + C; 2m + 1; 



2' ^ ' " " ' £ + 1 

x r(/3 + 2 + /I + C) 2*1 /? + 2 + £ + C; 2/2+ l; ^-y ) . (A..13) 

We are now ready to express 2 in its final form for n even, 
1 = 2 + - (-l) n X- 

m ^ _ _ c My _, +n 

x [sinvr(a + /i + C) B'(a, //, g% £, C) M, 9, T, C) 

- sinvr(a - M + C) B'(a, — g*, £, C) 9, t , ] , (A.14) 

where we have introduced the dimensionless parameter r q = Am 2 /\q\ 2 and the conformal 
blocks 

B ^^)-(jTlJ (^J 2 r^ + i) 

2 Fi M + 2 + a + C; 2 M +1; j^-jj ■ (A.15) 

We can make the replacement 

B'(a, n, q*,£, C) - B"(a, M, ?*, C), (A.16) 

where 

2Fi ^ + 2+ a+c / -i-„4 , ) , +1 ,_/') i (A17) 

by using the identities 9.134:1 and 9.131:1 of [26] and assuming n is even, considering that 
the conformal block will be multiplied with another one. This displays explicitly that the 
£ <-> — £ symmetry depends on a and /?; there will be an overall factor in the product of 
conformal blocks, 

(f y +fi (fr) / * f " = (sign er_/3+n ' (A - 18) 

when analytically continuing into the physical region — 1 < £ < 1. This factor is pulled out 
of the conformal blocks in equation ( 4.22Q , appearing in the main body of the text. 
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The results ( A. 14 ) and ( A.17 ) show features of our earlier calculations of heavy [11] 
vector meson photoproduction for a = 0. For m — * 0, it gives the amplitudes for light 
vector meson production in the massless case [27|. In Appendix |B| we will show the result 
in the massless quark limit, and how to obtain it from the present massive quark result. 

Note that when evaluating these expressions numerically it is necessary to carefully 
consider the branch cut structure of the conformal blocks, which are not single- valued 
individually. It is then useful to further transform Eqs. (A. 15) or (A.17) by e.g. using 



formula 9.132:2 of [26], isolating the branch cuts outside of the hypergeometric functions. 



B. The massless quark limit 



In [27 1 the formulae corresponding to QA.14 ) and ( A.17| ) were computed in the approxima- 
tion of zero quark mass. This corresponds to using the small argument behaviour for the 
case a = 1, relevant for the calculation of chiral even amplitudes, 

mKi(m\p\) ~ l/\p\ (B.l) 

Then, one obtains a simpler result J = I\ m= o not involving the line integral over £: 



J 



\ 4 

(-i) a - p+n [sin7r(a + M-l/2)C(a,/i,g*,OCG9,At,9,r 
(-1)" sin n(a - p - 1/2) C(a, —p, q*,£) C((3, -p, q,C) ] , 



(B.2) 



where the C(a, p, q* , £) conformal blocks are given by 5 B"(a,p,q*,£,( = — 1/2) 



C(a,p,q*,{,) 



f) 



Ct+fl 



(C 2 



-(M+3/2+o)/2 



T(p + 3/2 + a) 



, > + 3/2 + q p-l/2-a _ _ 
2-^1 I ^ > n ! M + 1 ) 



rv + i) 
i 



i-e 



(B.3) 



The calculation proceeds in a similar way to the preceding one and we do not show the 
details here. It should, however, be possible to recover the massless result from the massive 
one by taking the limit m — > 0. This can be done by taking the leading pole approximation 
for the C integral in the limit r q = 4m 2 /\q\ 2 — > 0. The product of Gamma functions in the 
integrand becomes r(l/2 — C)T( — 1/2 — C) and there is a simple pole at £ = —1/2 and double 
poles at C = 1/2, 3/2, .... The double poles, however, vanish in the limit r q — ► 0, and using 
the residue theorem for the pole at £ = — 1/2 reproduces the massless result ( |B . 2|) . Note 
that formula (|B.2j ) is a generalization of an integral calculated in [2C] for a = /3, and we 
have checked that our result agrees. 

The perturbative chiral-odd photon wave function corresponds to a = and is propor- 
tional to mKQ{m\p\) — > — mln(m|p|), thus it vanishes in the massless limit. 

We stress that the above formulae (i.e. in the massless limit) are only reliable away 
from the regions of u close to or 1. This is because the substitution |B.1| breaks down 



5 The substitution £ = —1/2 is made to keep the definitions of a and (3 the same, because there is a 
difference by a factor of \p\ between the massless and massive integrand. 
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as p — > oo. In effect the expansion parameter is more like m 2 /(u(l — u)\q\ 2 ) and this 
is not small at the endpoints. In particular, a careful treatment reveals that there are 
no endpoint singularities in the full amplitude in the massless limit, for any individual 
conformal spin. In the z — » (two-gluon) limit, the endpoint singularities shown by || 
in the M|T n amplitude are recovered by a non-convergent sum over all conformal spins, 
analogous to that demonstrated in [18] for qq — > qq scattering. 
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